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Abstract 

In this paper we apply a variant of Heisenberg's quantization method for 
strongly interacting, non-linear fields, to solutions of the classical Yang-Mills 
field equations which have bad asymptotic behavior. After quantization we 
find that the bad features (i.e. divergent fields and energy densities) of these 
solutions are moderated. From these results we argue that in general the 
n-point Green's functions for Yang-Mills theories can have non-perturbative 
pieces which can not be represented as the sum of Feynman diagrams. A 
formalism for dealing with these non-Feynman pieces via nonassociative field 
operators is suggested. These methods may also find some application in 
dealing with high-T c superconductors. 
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I. INTRODUCTION 



The quantization of strongly coupled, nonlinear fields is a difficult and still open prob- 
lem. The success of QED and the Standard Model of electroweak interactions is largely due 
to the small size of the coupling constants, which allows for the perturbative expansion of 
physical quantities in terms of a power series in the coupling constants. This Feynman dia- 
grammatic technique is the standard method for dealing with weakly coupled theories. For 
strongly interacting field theories such as quantum chromodynamics (QCD) this procedure 
can not be applied in energy regimes where the coupling constant is large. The theoretical 
demonstration of confinement and the determination of the mass spectrum of light mesons 
and baryons are two examples where the perturbative techniques of quantum field theory 
are not good tools. These low energy QCD phenomenon must be dealt with as completely 
nonperturbative, quantum effects. 

A mathematical way of stating this basic difficulty in QCD is that one needs some 
nonperturbative method for calculating the n-point Green functions of the theory. Usually 
these problems in strongly coupled field theories are dealt with using numerical simulations, 
but the computing power of modern computers does not yet permit, in general, the same kind 
of precise comparison between theory and experiment as is found in QED or the electroweak 
theory. Even if the state of the art in lattice gauge theory calculations increases to the same 
level of precision, it would still be desirable to have approximate, analytical methods for 
investigating these non-perturbative aspects of QCD or any strongly coupled field theory. 
Here we apply a nonperturbative method for calculating the n-point Green's functions of a 
theory which was originally proposed in the 50's by Heisenberg (see Refs. [0, @ and ||) to 
deal with the strongly coupled, non-linear Dirac equation. In this method one writes down 
an infinite system of equations which connects together all the n-point Green's functions of 
the theory. Then by introducing some physically reasonable approximations one reduces this 
infinite system of equations into a finite system. Here we apply this quantization procedure 
to several classical solutions in SU(2) and SU(3) gauge theories. Although these classical 
solutions possess some interesting properties they also have bad features which are eliminated 
or reduced after the Heisenberg quantization is applied. 



II. CLASSICAL SINGULAR SOLUTIONS IN SU(2) AND SU(3) GAUGE 

THEORIES 

In this section we examine several classical solutions for SU(2) and SU(3) Yang-Mills 
gauge theories. These solutions all have similar, qualitative asymptotic behavior: one part 
of the classical gauge potential diverges as a power law and another part oscillates strongly. 
In this section we follow the notation and conventions of Ref . || • 

A. The SU(2) gauge "string" 

Consider the following cylindrical ansatz for SU(2) Yang-Mills theory 

A] = f( P ), (la) 
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A\ = v(p), (lb) 

4 = p w (p)i ( lc ) 

z, p, 4> are the standard cylindrical coordinates; a = 1, 2, 3 are SU(2) color indices. Inserting 
this ansatz into the Yang - Mills field equations 

1 -_d, (V^F a » u ) + f abc F b ^Al = 0, (2) 



yields 



f"+ f - = f (*" + ^) , (3a) 
v" + - = v i-f 2 + w 2 ) , (3b) 



u; + r = w (-/ +f 1 , (3c) 

p p z K ' 



In the case where w = Eqs. (0) reduces to 



/" + - p = fv 2 , (4a) 
v" + V - = -vf 2 . (4b) 



The asymptotic behavior of the ansatz functions /, v and the energy density £ are 

„2 



/ 



cos (2x 2 + 2<pi 
x H 



(5a) 



16x 3 

t, - y^ Sin(x2 + 0l) , (5b) 

x 

£ oc f 2 + v' 2 + fv 2 « const, (5c) 



where x = p/po is a. dimensionless radius, and po,0i are constants. A numerical study of 
Eqs. (H) was carried out in Ref. |4j], and it was found that A] = f(p) was a confining, 
power law potential and A 2 = v(p) was a strongly oscillating potential . Depending on the 
relationship between the initial conditions v(0) and /(0) the energy density near p = was 
either a hollow (i.e. an energy density less than the asymptotic value) or a hump (i.e. an 
energy density greater than the asymptotic value). On account of this and the cylindrical 
symmetry of this solution we called this the "string" solution. 

B. SU(3) "bunker" solution 

Next we considered the following SU(3) ansatz for the gauge fields 0: 

A ° = - *V + + ^ {K 3 + A %) ^>(r), (6a) 

At = (Afc - X%) ^ (/(r) - 1) + A? fc (c^x* + e ilk x j ) ^v(r), (6b) 
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here A a are the Gell - Mann matrices; a = 1,2,. ..,8 is a color index; the Latin indices 
i,j,k,l = 1,2,3 are space indices; i 2 = —1; r,9,4> are spherically coordinates. Substituting 
Eqs. (^) into the Yang - Mills equations (0) gives the following system of equations for 
f(r), v{r), w{r) and <p(r) 



r 2 f = f-f + 7fv' + 2vw<p -flw* + <p 



2 // 

r v 



r 2 w" 



v + 7vf 2 + 2fw(p - v (w 2 + y? 2 ) 



6w (/ 2 + 



Ufvcp, 



V" = 2<p (/ 2 + v 2 ) - Afvw. 



(7a) 
(7b) 
(7c) 
(7d) 



This set of equations is difficult to solve even numerically so we investigated simplified cases. 
First we examined the / = tp — case, under which Eqs. (0) became 



r 2 v" 
r 2 w" 



v 3 — v — v w 2 , 
6wv 2 . 



In the asymptotic limit r — > oo the form of the solution to Eqs. 
similar to the "string" solution: 



3A 



v ~ 
w ~ 

2 _ 



A sin (x a + O ) , 

a — 1 cos {2x a + 



± 



ax" + 
I). 



x u 



(8a) 
(8b) 

approaches a form 

(9a) 
(9b) 
(9c) 



where x = r/r is a dimensionless radius and r ,<j) , and A are constants. A numerical 
investigation of Eqs. ([]) was done in Ref. ||. Since t>(r) was strongly oscillating this 
resulted in the space part of the gauge field of Eq. fl6b|) being strongly oscillating. The 
ansatz function w(r) increased as some power of x as x — ► oo, which indicated a classically 
confining potential. 

The "magnetic" and "electric" fields associated with this solution can be found from A a , 
and have the following behavior 



Et oc 



E a r cx 



v 2 -l 
rw' — w 



m oc v \ 



vw 



HS oc t»', 



(10a) 
(10b) 



for E r a , if e a , and ^ the color index is a = 1, 3, 4, 6, 8 and for H?, E a e and E% a = 2, 5, 7. The 
asymptotic behavior of ff$, ifg and £?^, is dominated by the strongly oscillating function 
v(r). Later we will show that the quantum corrections to this solution tend to smooth out 
these strongly oscillating fields. ^From Eqs. QTO ) and the asymptotic form (Q) the radial 
components of the "magnetic" and "electric" have the following asymptotic behavior 



' »,2— a 



(11) 
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The radial "electric" field falls off slower than 1/r 2 (since a > 1) indicating the presence of a 
confining potential. The 1/r 2 fall off of H® indicates that this solution carries a "magnetic" 
charge. This was also true for the solutions discussed in Refs. [|J 0. 

This solution of the classical SU(3) field equations exhibits fields which lead to a classical 
confining behavior, and which have similarities with certain phenomenological potential 
models of confinement. The most significant draw back of these solutions is that they have 
infinite field energy. In the case of the "bunker" solution one finds that the asymptotic form 
of the energy density goes as 




2 a 2 (q-l)(3a-l) 



As a > 1 this energy density will yield an infinite field energy when integrated over all space. 
This can be compared with the finite field energy monopole and dyon solution || . However 
it has been demonstrated || that the finite energy monopole solutions can not trap a test 
particle while the infinite energy solutions can. 



C. The classical f,(p^0 case 
Finally we examine the v — w — case: 

r 2 f" = f~f- ftp 2 , (13a) 
r\" = 2 V f 2 (13b) 

These equations have well known, nonsingular monopole solutions, but these monopole 
solutions require special boundary conditions. If one does not impose special conditions 
then a more general solution is 

/ « A sin (x a + O ) , (14a) 
a - 1 cos (2x a + 20o 



(fi « ± 



ax" + 



(14b) 



4 X" 

A 2 = a(a-l). (14c) 

which is the same as the asymptotic form of the "bunker" solution, and thus shares its bad 
long distance behavior. A numerical investigation of the solutions to Eq. ( p~3| ) also indicates 
a similarity to the "bunker" solution. The "magnetic" and "electric" fields associated with 
this solution are also similar to those of the "bunker" solution. 

Ha r^ f -^i H;<xf, HS<xf, (15a) 
£?oc^P^, ^oc^, ^oc^, (15b) 

here the color indices a = 2, 5, 7 indicate that this is really an SU (2) gauge potential embed- 
ded within SU(3). We now review the Heisenberg quantization method in order to apply it 
to these three solutions. 
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III. HEISENBERG QUANTIZATION OF STRONGLY INTERACTING FIELDS 



Heisenberg's basic idea is that the n-point Green's functions can be found from some 
infinite set of coupled differential equations, which are derived from the field equations for 
the field operators. As an example we show how this method of quantization works for a 
spinor field with a nonlinear self interaction ]l| @ 0- The basic equation for the spinor 
field is : 

Yd^(x) - l 2 $$(x)$(x)$(x))] = (16) 

where 7^ are Dirac matrices; ^(x), ip(x) are the operators of the spinor field and its adjoint 

respectively; ^[^(ip'ip)] = ipfyip) or ?/>7 5 (?/>7 5 ?/>) or ^7^ (^7^) or ?/>7 At 7 5 ('?/>7 At 7 5 '?/>). The con- 
stant I has units of length, and sets the scale for the strength of the interaction. Heisenberg 
emphasized that the 2-point Green's function, G 2 (x 2 , ^1), in this theory differs strongly from 
the propagator in a linear theory in its behavior on the light cone : in the nonlinear theory 
G 2 (x 2 ,Xi) oscillates strongly on the light cone in contrast to the propagator of the linear 
theory which has a 5-like singularity. Heisenberg defined r functions 

r{x l x 2 ...\y l y 2 ...) = (0|T[^(x 1 )^(a: 2 )...^(i/ 1 )^(y 2 )...]|$) (17) 

where T is the time ordering operator; 1$) is a state for the system described by Eq. (|16|). 
Eq. (|TTD establishes a one-to-one correspondence between the system state, |$), and the set 
of functions r. This state can be defined using the infinite function set of Eq. (pj]). Using 
equation (|T6| ) and ( |T7j ) we obtain the following infinite system of equations for the r's 

d 

l ~ 2 rfr)^-lT T ( X l--- X n\yi---yn) = Q[r(x 1 ...X n X r ^...ynVr)} + 

(r) 

5{x r - yi)T(x 1 ...x r -ix r+1 ...x n \y 2 ...yr-iyr+i—yn) + 

K x r ~ y2)r{xi...x r ^ 1 x r+1 ...x n \yiy 2 ...y r ^iy r+1 ...y n ) + ... (18) 

Eq. ( |TBD represents one of an infinite set of coupled equations which relate various orders 
(given by the index n) of the r functions to one another. To make some head way to- 
ward solving this infinite set of equations Heisenberg employed the Tamm-Dankoff method 
whereby he only considered r functions up to a certain order. This effectively turned the 
infinite set of coupled equations into a finite set of coupled equations. 

Heisenberg used the procedure sketched above to study the Dirac equation with a nonlin- 
ear coupling. Here we apply this procedure to nonlinear, bosonic field theories such as QCD 
in the low energy limit. In particular by applying this method to the infinite energy, solutions 
just discussed, and making certain assumptions analogous to the Tamm-Dankoff cut-off, the 
unphysical asymptotic behavior of these classical Yang-Mills solutions is "smoothed" out. 
It is also possible that this quantization method could be applied to the recently proposed 
|T0| strongly interacting phonon theory of high-T c superconductors. 

IV. QUANTIZATION OF SINGULAR CLASSICAL SOLUTIONS 

Before applying the above quantization method we make the following simplifying as- 
sumptions : 
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1. The degrees of freedom relevant for studying the singular solutions (both classically 
and quantum mechanically) are given entirely by the two ansatz functions /, v for 
Eqs. OH) or, v,w for Eqs. (ffl) or f,ip for Eqs. (Jl3|). No other degrees of freedom arise 
through the quantization process. 

2. ^From Eqs. ([|) and (|TJ]) one function is a smoothly varying function for large x, 
while another function is strongly oscillating. We take the smoothly varying function 
to be an almost classical degree of freedom while the oscillating function is treated as 
a fully quantum mechanical degree of freedom. Naively one might think that in this 
way only the behavior of second function would change while first function stayed the 
same. However since both functions are interrelated due to the nonlinear nature of 
the field equations we find that both functions are modified. 

A. Pedagogical example : Weak anharmonic oscillator 

As a pedagogical example we repeat some details of Heisenberg's quantization method 
as applied to the weak anharmonic oscillator . This section is also intended as a reference 
for when we apply this same procedure to the various classical Yang-Mills solutions. The 
Hamiltonian equations for the operators of the weak anharmonic oscillator are 

q = ~, (19a) 
m 

I- = - u 2 q - \f (19b) 
m 

here q and p are the operators of the canonical conjugate variables for the anharmonic 
oscillator. The r functions are 

T (k\ n _ k) e ^ = (OlTqihMh) ■ ■■q{t k )p{t k+1 )p{t k+2 ) ■ • -p(t„)|$)| tl=t2= ... =tn= , (20) 
The Hamiltonian operator equations yield the following system: 

k 

iujT{k\l) = —r(k - l\l + 1) - mujlh{k + l\l - 1) - 



rn 



. 2 l(l-l)(l-2) 



mXlr(k + 3\l - 1) + mXh 2 ^ ^ J -r(k + l\l - 3) (21) 



To a first approximation this reduces to : 

WllO) = — r(0|l), (22a) 
m 

— r(0|l) = -u#r(l|0) - Ar(3|0) (22b) 
m 

This is not a closed system, but by making the following approximation: 

T(3|0) = li T(1|0) ' (23) 
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one can solve Eqs. (|22|) to obtain: 

3 h\ 



u = ±u Q [l + - 3 . (24) 



Which coincides with the quantum mechanical solution as an expansion of order A : uox,o = 
{Ex - E )/h 

B. The quantized SU(2) gauge "string" 

To apply Heisenberg's quantization scheme to the Yang-Mills system of Eq. (ffl) we 



replace the ansatz functions by operators f(p),v(p) as in [O 



f" + - = fv 2 (25a) 

x 

v" + - = -vf 2 (25b) 

x 

the primes denote derivatives with respect to the dimensionless radius x. Taking into account 
assumption (2) we let / — > / become just a classical function again, and replace v 2 in Eq. 
( p5a|) by its expectation value 

f" + - = f(v 2 ) (26a) 
x 

v" + - = -vf 2 (26b) 

x 



Now if we took the expectation value of Eq. (|26b| ) and ignored the coupling to / on the 
right hand side we would have an equation for determining (v) = (<&|?)|$). However the two 
nonlinear terms on the right hand side of Eqs. (|26a| - |26b| ) show that a new object, (v 2 ) 



enters the picture so that Eqs. (|26a| - 126 b| ) are not closed. To obtain an equation for (v 2 ) 



we act on v (x) with the operator + -^J giving 

(v 2 )" + ^(v 2 )' = -2v 2 f 2 + 2(v') 2 (27) 
Taking the expectation value of this equation gives the desired equation for (v 2 ) 

( v y + I(^y = - 2 (v 2 )f 2 + 2((v') 2 } (28) 

x 

Again this equation is not closed due to the ((v') 2 ) term. We could again try to find an 
equation for ((v') 2 ) by the same procedure we employed for (v 2 ). This equation would 
also not be closed. Continuing in this way we would find an infinite set of equations. In 
order to have some hope of handling this problem we need to make some approximation 
to cut this process off. For the anharmonic oscillator Heisenberg solved this by assuming 
that t(xx ■ ■ - x\ pp - ■ - p) = (x(t)x(t) ■ ■ ■ x{t) p(t)p(t) ■ ■ - pit)) w for large m or n. We try 
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two different approximations for the ((V) 2 ) term and show that both yield similar large x 
behavior that fixes the infinite field energy problem of the classical solution. First we assume 
that ((V) 2 ) ~ ±(v 2 )'. Since ({v') 2 ) is positive definite one picks the ± sign so that the right 
hand side of this assumption is also positive definite. Under this assumption the equations 
become 

(v 2 )"+(±T2)(vy = -2(v 2 )f (29a) 

f" + -f' = f(v 2 ) (29b) 

x 

The approximate solution of Eqs. (|29a| - |29b|) is of the form 

{y2) ^ vg exp(- 7 x) (3Qa) 



expj-'yx) 

J « too + Jo 7= (30b) 



where 

/o7 2 = /oo^o 2 7 2 ± 2 7 = -2/1 (31) 
The second relationship can be written (using the first relationship) as 

7 = T/oo (foe + J- J (32) 



Although (v' 2 ) ~ +(u 2 )' leads to unphysical exponentially growing solutions, the assumption 
(v' 2 ) pa —(v 2 )' leads to exponentially decaying solutions. Thus under this latter assumption 
we find that the quantum mechanical treatment of this nonlinear system modifies the bad 
features of the classical solution. The asymptotic behavior of v goes from being strongly 
oscillating (see Eq. (|5T5D) to decaying exponentially, while the asymptotic behavior of / 
goes from being linearly increasing (see Eq. fl5aD ) to also decaying exponentially. In some 
intermediate region the two solutions should match up with one another. If the asymptotic 
form for these ansatz functions are used in the energy density, £, of Eq. ( j5cD we find that 
the field energy is now finite. To calculate £ we would replace the classical terms v' 2 , f 2 v 2 
by the appropriate quantum operator and take the expectation value. The (v' 2 ) term would 
be handled according to the assumption we used for closing the equations. 



C. Quantized "bunker" solution 

Although the classical behavior of "bunker" solution is interesting due to its similarity 
with certain phenomenological confining potentials, the infinite field energy discussed at the 



end of section |IIB| strongly argues against the physical importance of this solution. One 
possible escape from this conclusion is if quantum effects weakened or removed this bad long 
distance behavior. 



Again in Eq. (g) we replace the ansatz functions by operators w(x),v(x) [jO 
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2 -// "3 - - - 2 /on \ 

XV = V — V — VW {66a) 

x 2 w" = 6wv 2 (33b) 

here the prime denotes derivatives with respect to the dimensionless radius x. Taking into 
account assumption (2) we let w — > w become just a classical function again, and replace v 2 
in Eq. ( |33b| ) by its expectation value to arrive at 

x 2 v" = v 3 — v — vw 2 (34a) 

x 2 w " = 6w(v 2 ) (34b) 



If we took the expectation value of Eq. ( 34a ) we would almost have a closed system of 



differential equations relating w and (v). The (v 2 ) term from Eq. ( ^4b| ) and the (v 3 ) term 
from Eq. ( |34a| ) prevent the equations from being closed. Applying the operation x 2 d 2 /dx 2 



to the operator v 2 and using Eq. fl34a|) yields 

x 2 (y 2 )" = 2v 2 (y 2 -l-w 2 ) + 2x 2 (v') 2 (35) 

If we took the expectation of the above equation with respect to fluctuations in the ansatz 
function operator v 2 , and combined this with Eq. ( |34b|) we would almost have a closed 
system for determining w and v 2 except for the ({v') 2 ) term on the right hand side of 
Eq. fl3"5]). Continuing in this way one could obtain an infinite set of equations for various 
powers of the ansatz function operator (i.e. v n ). As in the previous section we make some 
assumption that effectively cuts off the system of equations at some finite order. We assume 
that the mean square of the <fi component of the magnetic field = £>' is small [] 

[(AH*)' 2 ) = ((#*- (^)) 2 ) = (H}) - (h^) 2 « (36) 
that implies 

H})^{H,) 2 or ((v'f)*mT (37) 

Then by taking the expectation of Eq. ( pq) we arrive at a closed system of equations from 
Eqs. O (|4§ (0) 

x 2 (v 2 )" = 2(v 2 ) 2 - 2(v 2 )w 2 - 2(v 2 ) + 2x 2 ((v)') 2 (38a) 
x 2 w " = 6w(v 2 ) (38b) 
x 2 (v}" = (v) 3 - (v) -w 2 (v) (38c) 

here we have further assumed that (t> 4 ) w (v 2 ) 2 and (v 3 ) ~ (v) 3 . Note that the present 
situation is slightly more complex than the quantization of the string solution of the previous 
section since now we have three coupled equations rather than two. It is straightforward to 



show that in the limit x — > oo the closed system given by Eqs. ( |38a|) - (38c) is solved by 



1 in other words the fluctuations of the Ha component of the magnetic field are weak. 
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(v 2 ) « 1 + ^ (39a) 

to » (39b) 

x l 

(v) w ±1 + - (39c) 
x 

Eqs. ( |39a| ) - ( |39cj ) provide information about the behavior of the "classical" ansatz function 
w, and the "quantum" ansatz function, v, via (v) and (v 2 ). The main point is that after 
applying the Heisenberg-like quantization procedure to the classical singular solution, the 
infinite increase of the ansatz function w, has changed to an acceptable asymptotic behavior 
(i.e. one that leads to a finite field energy). By replacing v 2 , w, {y') 2 in Eq. (|i~2"D with 
(v 2 ) ,w, ((v)') 2 - from Eqs. (|39a| ) - ( |39c| ) - we find that the field energy density of the 
quantized "bunker" solution takes the form 

£ oc I (40) 

in the limit in which quantum fluctuations become important (i.e. for non-Abelian theories 
which exhibit asymptotic freedom this means in the low energy or r — > oo range) the energy 
density goes from the form given in Eq. (|TJ) to that given in Eq. ([40]) . This gives a finite 
field energy. In the high energy or short distance regime we assume that the fields approach 
the classical configuration of Ref. [[| due to asymptotic freedom. This classical configuration 
is well behaved at r = 0, but has an infinite field energy due to its divergence as r — ■> oo. 



D. The quantized f,ip^0 case 

Finally, we employ the above quantization technique to the singular solutions of Eqs. 
(|13|) for the v = w — case. Again we replace the ansatz functions by operators f(x), y3(x) 

x 2 f" = f 3 -f- /> 2 (41a) 
x 2 (p" = 20f 2 (41b) 

using assumptions similar to the previous sections we have the following equations for the 
quantum function /, and the classical function ip: 

x 2 f" = f 3 -f- f& 2 ) (42a) 
x 2 V " = 2v(f 2 ) (42b) 

Using the same assumptions (including the assumption about the weak fluctuation) and 
procedure as in the previous sub-section we can turn Eqs. ([42a ) into a closed system of 



equations for the expectation values of various powers of the ansatz function operators 

x 2 (f 2 )" = 2(/ 2 ) 2 - 2(/ V - 2(f 2 ) + 2x 2 ((f)') 2 (43a) 
*V = 2^(f 2 ) (43b) 
x 2 (f)"=(f) 3 -(f)-V 2 (f) (43c) 
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This system has the following asymptotic solution 



(f) « 1 - ~ (44a) 
x 

(/) « ±1 - ~, (44b) 

x 

b , v 

y? « - (44c) 



The field energy density for this solution takes the form 



£ oc i (45) 
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V. OPERATOR QUANTIZATION FOR STRONGLY INTERACTING FIELDS 

In the previous sections we saw that strongly nonlinear fields {e.g. non-Abelian 
gauge fields) can have nonlocal objects (the classical SU(2) flux tube solution and the 
SU(2),SU(3) "bunker" solutions) which have good long distance behavior after an appli- 
cation of the Heisenberg quantization procedure. This result is markedly different from the 
situation for linear quantum field theories with no or small self-interaction. In these cases 
the quantum fields can not form static field configurations such as monopoles, flux tubes, 
or "bunker" solutions discussed above. In general one needs strongly interacting, nonlinear 
theories in order for such objects to exist. 

All of this may be an indication that the n-point Green's functions in non-Abelian 
gauge theories can have some pieces which are not simply the sum of Feynman 
diagrams. Physically, this would imply that some physical processes in non-Abelian gauge 
theories can not be explained solely in terms of the perturbative fluctuations of field quanta. 
Heisenberg may have had this situation in mind when he indicated ]lj that the divergence 
behavior in a theory with strong interactions will differ from that found in a perturbatively, 
renormalizable theory like quantum electrodynamics. In this section we want to first present 
some formal and then some heuristic arguments which support this idea of non-perturbative, 
non-Feynman contributions to the Green's functions. 



A. Quantum field theory with nonassociative field operators 

In this sub-section we consider field operators which are nonassociative. There are two 
reasons for doing this. First, in the Heisenberg quantization scheme one problem (which was 
dealt with by an approximate Tamm-Dankoff cutoff) was the appearance of an infinite system 
of differential equations with interrelated Green's functions. Nonassociative field operators 
can lead to a quantum theory where the Green's functions are functionally independent 
of one another, which would help this problem. Second, a nonassociative quantum theory 
can lead to amplitudes and processes which cannot be represented simply in a Feynman 
diagrammatic way by the exchange of virtual quanta. This supports the general picture 
given in the introduction to this section. In the following we will use the notations and 



conventions of Ref. 13 



12 



In the standard, noncommutative quantum field theories it is impossible to define an 
n-point Green's function that is not some polylinear combination of propagators and inter- 
action vertices. This general feature of the Green's functions is altered if the field operators 
are nonassociative. If one considers the time-ordering of a product of n field operators then 
the time-ordering will result in a shifting of the positions of the individual field operators. 
In a noncommutative theory the shifting in the positions of the field operators usually leads 
to the appearance of a commutator. The time-ordering can also lead to a change in the 
ordering of the brackets which group field operators together. In a nonassociative theory 
such a change in the ordering of the brackets leads to the appearance of an associator, just 
as in a noncommutative algebra the displacement of an operator to the left or right gives 
rise to a commutator. A good feature of nonassociative quantum field theories is that one 
can obtain n-point Green's functions that are independent of each other. 

We denote the product of n field operators in which the brackets are arranged in accor- 
dance with some rule P by M n (P). For example, 

M 3 (P) = (46) 

where ip x = (p(x) is the field operator at the point x; in this case, the rule P is that all 
opening brackets are at the extreme left hand position. 

We call this ordering of the brackets for any number of operators, ( i. e. with all opening 
brackets in the extreme left hand position) the normal position ordering of the brackets, and 
denote it by colons 

:M n (P):= (^•••(y3 1 (vV--)y^, (47) 

On the left there are n opening brackets in succession. 

Next, we consider a nonassociative quantum field theory with the following axioms. 

1. Two monomials M n (Pi) and M n (P 2 ) that differ only in the placement rule of the 
brackets differ from each other by a numerical function A(xi,x 2 , ■ ■ ■ , x n ; Pi, P2), which we 
call the n-point associator: 

M n (Px) - M n (P 2 ) = A(x u x 2 , ...,x n ; P u P 2 ). (48) 

2. The action of the product of two monomials on the quantum state vector \v) is defined 
as follows: 

(M^M^P^v) = M n (P 1 )(M fc (P 2 )|u)), (49a) 
H(M n (P0M fc (P 2 )) = ((u|M n (P 1 ))M fc (P 2 ). (49b) 

3. We have the usual commutation rules 

[<p x ,<P v ] = G(x,y), (50) 

where G(x, y) is the 2-point Green's function (propagator) which differs from the propagator 
of a linear theory. 

As an example of this formalism we examine the following simplified model, in which the 
action of any field operator on the vacuum state annihilates it 
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<p x \vac) = {vac\<p x = 0. (51) 

Now the 3-point associator is given by 

(0 x y )0 z = x (0 y z ) + A(x, y, z). (52) 

Using the rules introduced above for the action of the field operators on the vacuum state, 
we find that A(x, y, z) = 0. There are two 4-point associators given by 

[(<Px<Py)'Pz)<Pu = (0xfiy)(0z0u) + A^X, y\z, u) , (53a) 
0x{0y(0z<Puj) = {0x0y){0z0u) + A 2 (x, y \ Z, u) , (53b) 

Acting on each of Eqs. (|53a|) from the right and left with the vacuum state and using the con- 
dition (0), the 4-point associators A\ j2 can be expressed in terms of the vacuum expectation 
value of the monomial ((p x ip y )(ip z cp u ): 

A 1:2 (x,y\z,u) = A(x,y\z,u) = -{vac\(0 x y )(0 z u )\vac). (54) 

~ y^^z^ ~ m m (x M , y 11 , z^ and m m are 4-vectors). Using the commutation properties 
of the operators, the property of the vacuum state, and taking the points to be pairwise 
spacelike (x M ~ y M , z^ ~ u^) we can establish the following symmetry properties of the 
4-point associator: 

A(x, y\z, u) = A(y, x\z, u) = A(x, y\u, z). (55) 

The Green's function can now be defined as 

T(M n (P)) = M n (P) + G( Xl ,x 2 , . . . , x n , P), (56) 

where T denotes the ordinary time ordering of the operators in the monomial M n (P), and 
also normal ordering of the brackets in the monomial M n (P). Thus each n-point Green's 
function can be expressed in terms of a polylinear combination of n-point associators (m < n) 
and propagators. 



B. Heuristic argument for non-Feynman contributions to the Green's functions 

In carrying out our approximate application of Heisenberg's quantization method to 
certain classical, singular solutions, we examined only two degrees of freedom (/ and v , v 
and w or / and ip). In each of these pairs one function was taken as a quantum degree of 
freedom and the other remained a classical degree of freedom. To investigate this method 
further one could drop this assumption by treating both functions as quantum degrees of 
freedom. In this section we look at a broader issue of how to include the other frozen degrees 
of freedom. In other words the quantization process might lead to more degrees of freedom 
beyond the original classical degrees (e.g. f and v). 

Our main assumption in addressing this question is that for any n-point Green's function, 
for strongly interacting, nonlinear fields, there is a non-perturbative piece which can be 
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treated by the Heisenberg quantization method 0. In the classical region these degrees of 
freedom region can be associated with classical (possibly singular) field configurations. The 
others degrees of freedom are then handled using ordinary perturbative methods (i.e. using 
Feynman diagram techniques). 

Physically, this means that in a theory with strongly interacting fields one has: 

1. nonlocal static objects, formed by those quantum degrees of freedom which are con- 
nected with (possibly singular) solutions of the classical field equations. 

2. the quanta which fluctuate around these nonlocal objects. 

In terms of QCD this implies that there are degrees of freedom which play a dominate 
role in the Feynman path integral to a first approximation. These degrees of freedom are 
connected with classical (possibly singular) solutions of the Yang-Mills field equations. This 
is somewhat similar to the saddle-point method of evaluating the Feynman path integral. 
However, in the saddle-point calculation one approximates the path integral via the value 
of the action evaluated at a particular solution 



here <p is a field degree of freedom (in our case it is the non-Abelian potential A"), and (po is 
a classical, nonsingular solution of the field equations. For singular, classical solutions this 
approximation is not valid. Instead one can as a first approximation take 



here (p c i represents the ansatz functions used in the classical field equation. For the ansatz 
functions of Eq. (|7|) (f,v,w,ip) the integration would become Dip d — > D f ' DvDwDip. This 
leads to the following expression for the path integral 



here (Af,) , are the classical solutions in terms of the ansatz functions, and (Af) , are the 

x t 1 ' cl x dev 

perturbative, quantum deviations from the classical solution. The first integral is calculated 
by some nonperturbative method (lattice gauge calculation, Heisenberg method etc.). The 
second path integral can be handled by ordinary perturbative methods (Feynman diagrams). 

At this point one can ask how renormalization fits in with the above development where 
the Green's functions contain a non-perturbative piece. In the standard perturbative treat- 
ment of a quantum field theory one must introduce a renormalization prescription. The 
physical basis for this can to some extent be traced to the point-like nature of the quanta 
and interaction vertices. Mathematically the divergences come from momentum integrals 
over closed loops which have propagators (which are distributions in linear field theories) 
raised to some power. In the developments of the previous sections, where the n-point 



2 this piece may be connected with the nonassociative properties of fields operators 




(57) 




(58) 




(59) 
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Green's functions have a non-perturbative part due to the non-linear nature of the quan- 
tized fields, it is no longer necessarily true that the Green's functions need to be distributions. 
A generalization of the standard Feynman diagram technique would represent any complete 
n-point Green's functions as the sum of two pieces : one the standard, perturbatively calcu- 
lated field fluctuations, and the other would be the non-Feynman part (possibly related to 
nonsingular or weakly singular classical solutions). A similar situation arises in superstring 
theory where the building blocks for the scattering diagrams are nonlocal strings. In both 
the present case and in the case of string theory the divergences of the quantum theory are 
eliminated or decreased fj 



VI. DISCUSSION AND CONCLUSIONS 

The basic idea proposed here is that quantum field theories with strongly non- 
linear fields can have the non-Feynman pieces. Also the classical, possibly singular 
solutions can play an important role in the corresponding quantum field theory. After 
quantization these classical solutions can lead to nonlocal objects with good long distance 
behavior. These quantized nonlocal object are connected with the nonperturbative part of 
the Green's functions, and they are not the result of the perturbative quantum fluctuations 
of the field. 

This approach to the quantization of non-Abelian field theories is very similar to the 
attempts to quantize gravity on the basis of the Wheeler-de Wit equation. In the present 
case there is an important distinction : the operator of canonical momentum is not simply 
related to the derivative of a canonical coordinate (i.e. V ^ fr where q is a field operator). 
This follows from the fact that the propagator has an unusual behavior on the light cone. 
Reversing the analogy between the non-Abelian field equations and the Wheeler-de Wit 
equation one could postulate that the Wheeler-de Wit equation may be more complicated 
than is normally written, since the momentum operator in gravity is also not simply the 
ordinary derivative with respect to some corresponding field variable. 



In sections |IVQ and [IV D| we found two spherically symmetric nonlocal objects with good 
asymptotic behavior after quantization. It may be possible to forge a connection between 
these solutions and the QCD picture for the "dressing" of valence quarks by gluons and 



quark- ant iquark pairs. The authors of Ref. write: "One views dressing as surrounding 
the charged particle with a cloud of gauge fields. We . . . see that, . . . this cloud spreads out 
over the whole space, resulting in a highly nonlocal structure. This is not unexpected in 
QED, but in QCD we will see that there is a nonperturbative obstruction to the construction 
of this dressing." The Heisenberg quantization method that we have outlined might provide 
a nonperturbative description of the above "dressing" procedure. If this is correct then our 
results give structures which are "dressed" with gluons but not with quark- ant iquark pairs 



3 A related question would be to investigate whether there is a classical, nonlinear field theory 
which after quantization resulted in n-point Green's functions with non-perturbative pieces that 
had a structure similar to p-branes. If so then one might consider string theory as a quantum field 
theory arising from some nonlinear (possible nonlocal) field theory ! 
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since we did not consider quark fermion fields in this work. 



A. Connection with the Maximal Abelian gauge model 

The quantization method of non-Abelian gauge fields given in this work shares some 
common features with the Abelian Projection (AP) model of gluodynamics proposed by 
't Hooft in Ref. |15|| . The basic idea of the AP model is that in a non-Abelian gauge 
theory such as QCD there are special degrees of freedom which play an essential role in 
confinement. These degrees of freedom are connected with monopoles which appear in the 
maximal Abelian subgroup, H, of the full non-Abelian group G. In the AP model the initial 
SU(3) field equations are rewritten as an Abelian gauge theory with magnetic monopoles 
and some matter fields. Roughly speaking, the AP method implies that the off-diagonal 
gauge fields are excluded from consideration. More precisely one integrates over all the 
off-diagonal (G/H) gauge fields, leaving only the Abelian degrees of freedom |T5|. 



SU(3) gauge theory has monopole solutions, which at large distances look like U(l) 
Abelian monopoles (i.e. the solution exhibits 1/r 2 behavior asymptotically). Thus af- 
ter quantization these Abelian monopole solutions should play a more important role in 
comparison with the quantized solutions examined in sections [IV C| and |1V D| since these 
quantized solutions have a faster fall off. This is in agreement with the AP model. 



B. Connection with High-T c superconductivity 

The quantum theory of solids can be taken as a quantum field theory on some lattice. 
Thus one could ask if the above formalism has an application to condensed matter systems. 
For example can the quantized, classical solutions discussed above play a role in condensed 
matter ? In Ref. [17]] it was shown that a flux tube solution can appear in superconductors 
in the mixed state. In [JRJ a model of Cooper pairing was proposed in which the phonons ex- 
changed between Cooper electrons are confined in a "flux" tube similar to the chromoelectric 
flux tubes which are thought to occur between quarks and anti-quarks in QCD. 

One possible application of such a model []TBj] might be as a mechanism to explain the 
Cooper pairing in high-T c superconductors. One could postulate that the Lagrangian for 
high-T c materials had a strong, nonlinear potential term. This could result in a classi- 
cal, cylindrical solution (either singular or nonsingular) , which after quantization via the 
Heisenberg-like method could give rise to a phonon flux tube that stretched between two 
Cooper electrons, binding them to much higher temperatures than in ordinary BCS super- 
conductors. 
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